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ABSTRACT 

Aims. To calibrate the relationship between Ne302 (A'e302=log([Nem]/i3869/[Oii]/13727)) and oxygen abundances, we present a 
sample of -3000 Hn galaxies from the Sloan Digital Sky Survey (SDSS) data release four. They are associated with a sample from 
the literature intended to enlarge the oxygen abundance region. 

Methods. We calculated the electron temperatures (T e ) of 210 galaxies in the SDSS sample with the direct method, and T e of the 
other 2960 galaxies in SDSS sample calculated with an empirical method. Then, we use a linear least-square fitting to calibrate the 
Ne302 oxygen abundance indicator. 

Results. It is found that the Ne302 estimator follows a linear relation with 12 + log(0/H) that holds for the whole abundance range 
covered by the sample, from approximately 7.0 to 9.0. The best linear relationship between the Ne302 and the oxygen abundance is 
calibrated. The dispersion between oxygen abundance and Nei02 index in the metal rich galaxies may come partly from the moderate 
depletion of oxygen onto grains. The Ne302 method has the virtue of being single-valued and not affected by internal reddening. As a 
result, the Nei02 method can be a good metallicity indicator in the Hn regions and Hn galaxies, especially in high-redshift galaxies. 
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> 1. Introduction 

in : 

Q\ , Those hydrogen II galaxies with strong emission lines are im- 
O^l ■ portant probes for the formation and evolution of galaxies. Their 
""xt" ' spectra contain the much important information needed to deter- 
, mine the star formation rate, initial mass function, element abun- 
■ dance, etc. (Stasihska & Leitherer 1996; Kennicutt 1998; Contini 
£■ — | et al. 2002). The heavy element abundance is a key parameter 
C^) . for the formation and evolution of a galaxy. Oxygen is the im- 
IJ ' portant element that is easily and reliably determined since all 
. i-h ! the most important ionization stages can be observed. The oxy- 
• gen abundance from the measurement of electron temperature 
Sh ; from [O m],U4959,5007/[Oiii]/l4363 is one of the most reliable 
5^ , methods. But [Ora]/l4363 is usually weak in the low metallicity 
galaxies, and there are often large errors when measuring this 
line. In high metallicity galaxies, [Om]/l4363 are hardly even 
observable. 

Instead of the T e method, strong line methods, such as the 
^23 P EL N2 Q, or 03N2 methods, are used widely (Pagel 
etal. 1979; Kobulnicky, Kennicutt, & Pizagno 1999; Pilyugin 
etal. 2001; Chariot & Longhetti 2001; Denicolo etal. 2002; 
Pettini & Pagel 2004; Tremonti etal. 2004; Liang etal. 2006). 
The R23 and P methods suffer the double-valued problem, requir- 
ing some assumption or rough a priori knowledge of a galaxy's 
metallicity in order to locate it on the appropriate branch of the 
relation. The N2- and 03N2 methods are monotonic, but the 
reasons for this are not purely physical. It is partly due to the 
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1 tf 2 3=([Oii]/t3727+[Om]/t/t4959,5007)/He 

2 ,P=[Oni],U4959,5007/([O ii]/13727+[0 ra]/l/l4959,5007) 

3 A'2=log([Nii]A6583/Ho;) 

4 03A'2=log(([0 ra]/l5007/He)/([Nii]^6583/HQ')) 



N/O ratio increases on average with the increase in metallic- 
ity (Stasihska 2006; Shi etal. 2006). Besides, calibrations of the 
03N2 and N2 indices might be improper for interpreting the 
integrated spectra of galaxies because [Nn]/16583 and Ha may 
arise not only in bona fide H n regions, but also in the diffuse 
ionized medium. Stasihska (2006) has recently proposed A r3 03 
and 53(930 as new abundance indicators, which have the ad- 
vantage of being unaffected by the chemical evolution effects. 
The advantages are superior to previous N2 and 03N2 methods. 

In short, one ideal metallicity indicator has to be mono- 
tonic and be independent of the internal reddenning and chem- 
ical evolution. Nagao et al. (2006) find that the Ne302 index, 
i.e. log([Nem]^3869/[On]^3727), fulfills these requirements. 
They derive a relation of Ne302 vs. 12 + log(0/H) by using 
the Bayesian abundances provided by the MPA/JHU group0 for 
the metal rich galaxies in the SDSS (Tremonti etal. 2004) and 
the T e abundances based on an electron temperature-sensitive 
line, [Om]/l4363, for metal poor galaxies. The Bayesian metal- 
licity could be problematic in some aspects, such as improper 
subtraction of the underlying starlight or unsuitable treatment of 
nitrogen enrichment in the HII galaxy model. As a result, there 
is a significant systematic difference between Bayesian metal- 
licity and T e metallicity (Yin et al. 2007), so it is necessary to 
recalibrate the Ne302 index based on abundance derived from 
T e method, which is believed to be the most reliable method 
to determine oxygen abundance. Also, the dispersion between 
[Ne iii]/13869/[Oii]/13727 and oxygen abundance needs to be ex- 
plained. 



5 Ar3O3=[Arm]^7135/[Oiii]i5007 

6 S3O3=[Sm]^9069/[Om]^5007 

7 See http://www.mpa-garching.mpg.de/SDSS/ 
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In this paper, we recalibrate the Ne302 metallicity index 
based on a large sample of H n regions and H n galaxies where 
oxygen abundance is determined by the T e method. For low 
metallicity galaxies (Sample I, see Sect.2), we used a five-level 
statistical equilibrium model in the IRAF NEBULAR package 
(de Robertis etal. 1987; Shaw & Dufour 1995), which makes 
use of the latest collision strengths and radiative transition prob- 
abilities to determine the T e and n e . For high metallicity galaxies 
(Sample II, see Sect.2), an empirical relation of T e and strong 
spectral lines from Pilyugin (2001) was adopted for the electron 
temperature determination (their Eq. 11). We also study the phys- 
ical origin of the dispersion between [Neiii]/l3869/[On]/l3727 
and oxygen abundance. 

The Sloan Digital Sky Survey (SDSS) is the most ambitious 
imaging and spectroscopic survey to date, and it will eventually 
cover a quarter of the sky (York et al. 2000). The large area cov- 
erage and moderately deep survey limit of the SDSS makes it 
suitable for studying the evolution of galaxies. Because of its 
homogeneity, the SDSS provides a large sample of Hn galaxies 
where oxygen abundance can be calculated with the classic T e 
method. The sample can be used to calibrate some metallicity 
indicators. 

This paper is organized as follows. Based on an SDSS DR4 
starbursts spectral sample and a sample of H n regions or H n 
galaxies from the literature, we present a sample to use for our 
metallicity index calibration (Sectf2]i. In Sect. [3] we calibrate the 
Ne302 metallicity index. In Sect|4] we study the origin of dis- 
persion between oxygen abundance and Ne302 index and check 
the accuracy of our calibration. In Sect. [5] we summarize the 
calibration result and discuss the merits of then Ne302 indicator 
with respect to the other strong line indicators and then conclude 
this paper. 

2. Data sample 

The preferred method for determining oxygen abundances in 
galaxies is obtained from the ratio of auroral to nebular line in- 
tensities, such as [Oiii],L14959,5007/[Om],14363 (the so-called 
T e method). In this paper, the adopted metallicities of Hn re- 
gions and H n galaxies for calibration are determined from the 
T e method. 

We used H n galaxies from the Fourth Data Release (DR4) 
of the SDSS. After subtracting the underlying starlight using the 
method of Li et al. (2005) and Lu et al. (2006), we fit the emis- 
sion line using the method of Dong et al. (2005). We made the 
internal reddening correction for the flux of all the emission line, 
using the two strongest Balmer lines, Ha/H/3 and the effective 
absorption curve = Ty(/l/5500A)~ - 7 , which was introduced 
by Chariot & Fall (2000). Then, we made use of the spectral 
diagnostic diagrams from Kauffmann et al. (2003) to classify 
galaxies as starburst galaxies, active galactic nuclei (AGN), or 
unclassified. To reduce systematic and random errors from aper- 
ture effects, our galaxy samples are limited by the requirement 
that redshift z > 0.04 (Kewley et al. 2005). 

Within the primary starburst sample, two subsamples were 
selected from the SDSS-DR4 with the fluxes of spectral lines 
for all [Nera]^3869, [O n] ,13727, H/3, [Om],l4959, [Om],15007, 
Ha and [Nn]/16583 higher than 5 times the flux uncertainty. 
The difference of these two subsamples is that the first sub- 
sample (Sample I) was selected by the additional criterion from 
the [O in] ,14363 line, and the flux uncertainties for [Om] ,14363 
was higher than 5<x. In all, 210 galaxies were included in 
this subsample. The [Om]/l4363 is strongly dependent on the 
metallicity of galaxies and it becomes undetectable in high 



metallicity galaxies. Therefore, galaxies in Sample I are those 
with low metallicity. The electron temperature in the high- 
temperature zone (T e (0 III)) of Sample I galaxies was calculated 
from [Om]ii4959,5007/[Oin]^4363, using a five-level statis- 
tical equilibrium model in the IRAF NEBULAR package (de 
Robertis etal. 1987; Shaw & Dufour 1995). In the second sub- 
sample (Sample II), galaxies have weak or no [Ora]/14363 line, 
and 2960 galaxies were selected in this subsample containing 
generally metal rich galaxies. The r e (OIII) of the galaxies in 
Sample II was calculated with the empirical method following 
Pilyugin (2001) from the R23 and P parameter (their Eq.l 1). 

To estimate the temperature in the low-temperature zone 
r e (OII), the relation between T e (OlT) and r e (OIII) from 
Garnett (1992) are utilized : 

t e (Oll) = 0.7xf,(OIH) + 0.3, (1) 

where t e =T e /10 4 K. After calculation of T e , we used the expres- 
sions from Izotov et al. (2006) (their Eqs.3 and 5) to calculate the 
oxygen abundances. The average uncertainty from flux measure- 
ment in the computed 12 + log(0/H) values is typically 0.10 
dex in both samples. 

To widen the oxygen abundance range, we built a large 
database of published reddening-corrected line fluxes for Hn 
regions or Hn galaxies, besides the Hn galaxies from SDSS 
DR4. The data for Hn regions in spiral galaxies in the lit- 
erature with good [Nem]/13869, [On]/l3727 and reliable T e 
metallicity are taken from Garnett & Kennicutt (1994), Garnett 
etal. (1997), Bresolin etal. (2004), Bresolin etal. (2005) and 
van Zee etal. (1998). The data for Hn galaxies in the litera- 
ture with very good signal-to-noise ratios are taken from Guseva 
et al. (2003a,2003b,2003c), Izotov et al. (1994, 1997, 1998, 
2004a, 2004b, 2006), Lee et al. (2003a, 2003b), Lee et al. (2004), 
Melbourne et al. (2004), Papaderos et al. (2006), and Thuan et 
al. (1995). In total, 84 H n regions of spiral galaxies and 446 H 11 
galaxies in the literature are adopted in our sample. The average 
uncertainty from flux measurement in the 12 + log(0/H) values 
is typically within 0. 10 dex. 

3. Ne'iOl as a metallicity calibrator 

The upper panel of Fig.Q]shows the oxygen abundance as a func- 
tion of Ne302 for our samples. Apart from a few outliers in the 
SDSS DR4 sample, all the objects merge into a relatively tight, 
linear, and steep sequence, as expected. As expected from stel- 
lar nucleosynthesis (Izotov & Thuan 1999) and confirmed by 
observation of low metallicity HII galaxies in the SDSS DR3 
(Izotov et al. 2006), Ne and O are produced by the same stars, 
so the relationship between oxygen abundance and the Ne302 
index unlike 03N2 or N2 is not affected by the chemical evo- 
lution effects. The observed trend of Ne302 is due to the fact 
that the ionization parameter has a strong metallicity dependence 
(Nagao etal. 2006). The two emission lines, [Nem],l3869 and 
[On]/13727, have different degrees of ionization. Their ratio is 
strongly dependent on the ionization parameter and correlates 
with metallicity. 

To confirm this view, we plot the relationship in 
Fig. [2] between Ne302 index and the P parameter for 
Hn galaxies in SDSS DR4. The P parameter is defined 
as ([Oiii]ii4959,5007)/([On]^3727+ [Oiii],U4959,5007) by 
Pilyugin (2001), and it is a good representive of an ionization 
parameter. Figure [2] shows that the Ne302 index is positive re- 
lated to the P parameter, which supports the idea that the Ne302 
index is negatively correlated with the metallicity, since more 
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metal rich H n regions are excited by a softer radiation field and 
have a lower ionization parameter. We have checked the relation- 
ship further between electron temperature T e and Ne302 index 
in Fig. [3] and find the observed trend of Ne302 is partly due to 
the increase in T e as metallicity decreases, which leads to an in- 
crease in the Ne302 ratio. The comparison between Figs.[2]and 
[3] shows that, for the observed relation between 12 + log(0/H) 
and the Ne302 index, the contribution from the decrease in 
the ionization parameter P with the increase in metallicity is 
stronger than the decrease in T e with an increase in metallic- 
ity, since the Ne302 index is more sensitive to the P parameter 
than T e . 

From the upper panel of Fig. [T| we can define a new metal- 
licity calibration. The observed distribution of all the points in 
this figure is linear least-square-fitted by the following expres- 
sion given as the solid line in the lower panel of Fig. \T\ 

12 + log(0/H) = -1.171(±O.O88)xJVe302 

+7.063(+0.727)) (2) 

with a standard error of 0.302 dex. In this panel, we also give 
the mean of the data with bin=0.091 dex in 12 + log(0/H) (the 
stars) and the result of fitting the mean. We plot the calibration 
of Nagao et al. (2006) to compare with our result. In Sect.4.2, we 
discuss why we prefer the fit for all the points and do not use the 
fit of mean. 

4. Discussion 

4.1. The origin of the dispersion 

Although the T e metallicities are closely related to Ne302 in- 
dex, the dispersion of the relation is large, especially in the low 
metallicity region. Some of the scatter may come from the mod- 
erate depletion of oxygen onto grains in the metal rich galax- 
ies. To show this view clearly, we calculated the log(Ne/0) 
abundance ratios for Samples I and II using the method of 
Izotov et al. (2006) (their equations (7) and (19)). Then we plot- 
ted the difference between \og{0 1 H)j e -\og{0 1 H)^ e w2 an d the 
log(Ne/0) abundance ratios in the upper panel of Fig. [4] for 
metal poor galaxies and the lower panel for metal rich galaxies. 
To compare Sample I and metal-poor emission-line galaxies in 
Izotov et al. (2006), we show the sample of Izotov et al. (2006) 
in the upper panel of Fig. [4] 

As expected, Sample I is consistent with the sample of 
Izotov etal. (2006). There is a positive correlation between 
(\og{0 1 H) Tr )-(\og{0 1 H) NeW2 ) and the log(Ne/0) for metal-poor 
galaxies because more oxygen depletion would be present in the 
H ii regions with higher metallicity (Izotov et al. 2006 ). 

The lower panel of Fig. [4] shows more clearly that there is 
a clear correlation between log((9 / ' H)T t -\og(0 / 'H)^ e 302 an d the 
log(Ne/0) for the metal rich galaxies because more oxygen is 
locked in the dust grains in the more metal-rich Hn regions 
(Izotov et al. 2006 ). This correlation for Sample II can be given 
as a linear least-square fit: 

\og{OIH) Te - \og{0/H) Ne302 = 0.610(+0.024) + 

0.999(+0.037) x log(Ne/0), (3) 

with a standard error of 0.104 dex. The corrected calibration 
of Ne302 methods by considering this correlation will provide 
more accurate oxygen abundance, and the standard error of the 
corrected calibration in Fig.Q]will decrease to 0.198 dex, which 
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log(NelllX3869/0IIA3727) 

Fig.l. Upperpanel: 12 + log(0/H)vs.log(Ne302)forourdata 
sample. Dots are galaxies from SDSS DR4. Crosses are from the 
literature. Lower panel: The solid line represents Eq. (2). Stars 
denote the mean for each bin of the oxygen abundance. The dot- 
dashed lines denote the best-fit function for each metallicity bin. 
The dashed curves are the calibration of Nagao et al. (2006). 



is much less than the previous ones (0.302 dex). This corre- 
lation can be used to correct the derived oxygen abundances 
from Ne302 methods. However, to apply the Ne/O correction 
to 12+log(0/H) abundances, we have to measure electron tem- 
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Fig. 2. The relationship between Ne302 index and the P 
parameter for Hn galaxies in SDSS DR4. P=([Om]/l4959+ 
[O in] A5007)/( [O n] A3121 + [O ra]^4959+ [O m]/l5007) 
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Fig. 3. The relationship between Ne302 index and the electron 
temperature for Hn galaxies in SDSS DR4. The electron tem- 
perature is in units of 10 4 K. 



perature T e , H/3, [Om],t4959, [Om]^5007, besides [Nem]/13869 
and [On]/l3727, so it will not always be practical. 

4.2. Comparison with others 

We plot the calibration of Nagao et al. (2006) to compare with 
Eq. (2) in the lower panel of Fig. Q] Some samples of both 
ours and Nagao's are selected randomly from the literature. 
As a result, the calibration might be affected by selection ef- 
fects. To mitigate this effect, Nagao's calibration fits the mean 
in given metallicity bins. However, most of the SDSS data in 
our sample are concentrated in the metallicity range of 8.0< 
12 + log(0/H) <8.5. As a result, if we fit the mean, the calibra- 
tion will lose a lot of information in 8.0< 12 + log(0/H) <8.5, 
because thousands of SDSS data are represented by several dif- 
ferent mean values. As a result, the calibration fitting the mean 
will induce larger errors. In considering this, we used the calibra- 
tion fitting all data points in the lower panel of Fig. [T]to compare 
with Nagao's calibration. 

The lower panel of Fig. Q] shows that the oxygen abundance 
from Nagao's calibration is systematically higher than our cal- 
ibration at the given Ne302, except that at the low metallicity 
end, our calibration is consistent with Nagao's because both ours 
and Nagao's abundances at the low metallicity end are derived 
by the T e method. The difference between Nagao's calibration 
and ours in the high metallicity region could be caused by the 
oxygen abundances of most galaxies in the Nagao's sample be- 
ing calculated by the Bayesian methods, while the oxygen abun- 
dances of our sample are all derived with the T e methods. The 
origin of the difference between the Bayesian and T e metallic- 
ities have already been discussed by Yin et al. (2007). They 
find that, for almost half of the sample galaxies (227 among 
531 galaxies with T e measurements), Bayesian metallicities are 
overestimated by a factor of about 0.34 dex on average, which is 
consistent with our result. They propose that the overestimates of 
Bayesian metallicities may be related to the onset of secondary 
N enrichment in models. Another reason for the lower T e metal- 
licities than Bayesian metallicities is that the [Om]/l4363 emis- 
sion line is biased by the very hot H n regions in each galaxy; 
thus, the global average temperature might be overestimated to 
1000-3000K, which results in systematic underestimation of the 
oxygen abundance of 0.05-0.2 dex, as Nagao et al. (2006) pro- 
pose. 



i 

o 




-2 t , 

-1 .5 



.0 -O.b 0.0 

log(Ne/0) 



I 
en 

° 



2. -1 



-2 



•1.5 



-0.5 



0.0 



log(Ne/0) 



Fig. 4. Upper panel: Correlation between log(0/H)r e - 
log(0/ii)jVe3O2 an d the log(Ne/0) abundance ratios of Sample 
I (crosses) and the sample of Izotov et al. (2006) (dia- 
monds). log(0 / H)r r -log(0 / H)Neioi is the difference between 
T e metallicity determined from the expressions in Izotov 
etal. (2006) and Ne302 metallicity determined from Eq.(2). 
Lower panel: Correlation between log(0 / if)r e -log(0 / H);ve302 
and the \og(Ne/0) abundance ratios of Sample II (dots). 



In Fig|5] we compare the T e and Bayesian metallicities in 
our SDSS sample. It shows that the Bayesian metallicities of 
nearly half of the Sample II, where T e are derived from the em- 
pirical method in Pilyugin (2001), are ~0.3dex higher than T e 
metallicity. The relationship between the T e and Bayesian metal- 
licities behaves very much like the result in Yin et al. (2007) 
(their Fig. 2), which supports the empirical T e determination by 
Pilyugin (2001). This is consistent with T e derived by the classic 
temperature-sensitive line. Similar evidence can be found in Shi 
etal. (2006) (their Fig. 2) where T e metallicities from the empir- 
ical T e determination have nearly the same relation with other 
strong line metallicities as do those from T e derived from the 
classic temperature-sensitive line. 

Generally speaking, all the methods should result in the same 
abundance value for a given nebula, but this is not the case in 
practice. It is evident that this discrepancy is caused by the prob- 
lems both with models of HII regions and calibration. Why is the 
calibration of some types, such as the T e method in this paper, 
regarded as the most reliable way to derive the oxygen abun- 
dance? 

When comparing the numerical HII region models from 
Chariot & Longhetti (2001), which are the basis of the 
Bayesian abundances, with the numerical models of other au- 
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Fig. 5. Comparison of the T e metallcities and Bayesian metallci- 
ties. Crosses are Sample I. Small dots are Sample II. 



thors (Stasinska & Leitherer 1996; McGaugh 1991; CLOUDY 
by Ferland et al. 1998; or Kewley & Dopita 2002), one finds that 
there is significant disagreement between them because these 
models use different atomic data or different assumptions, and 
the stellar evolutionary synthesis code and photoionization code 
used in these models continuously improved. As a result, the 
present-day models cannot provide any uniform oxygen abun- 
dances. 

As for calibrations, the calibrations based on the measure- 
ments of real Hn regions is believed to be the most favorable 
way to derive the oxygen abundances, because it is calibrated by 
single H n region and will not be biased by other hot H n regions. 
That is why we use empirical T e metallicity rather than Bayesian 
metallicity in high metallicity regions to calibrate Ne302 index. 
There is a problem that the amount of H n regions with accurate 
measurements (a number of calibrating points) are not numerous 
enough. Therefore one has to use extrapolation, which can be a 
cause of uncertainty in the abundance determination. 

Because Ar and Ne are produced by the same star, as ex- 
pected by stellar nucleosynthesis (Izotov & Thuan 1999), it 
is instructive to compare the Ne302- and Ar303-methods 
(Stasinska 2006). We show this comparison for SDSS DR4 sam- 
ple in Fig. [6] 

Although having a large scatter and a clear tendency of the 
observed relation to slightly deflect from 1:1, there is agreement 
between the oxygen abundance from the Ne302- method and 
the Ar3(93-method for most galaxies. The dispersion between 
the r e -based and Ar303 methods is especially large in the mid- 
dle metallicity region (7.95 < 12 + log(0/H) < 8.2). It may 
be caused by the fact that, when calibrating Ar303 index, the 
oxygen abundance is calculated using the empirical P-method, 
not the r f -based method (Stasinska 2006), and the empirical P— 
method is problematic in the middle metallicity region, espe- 
cially 7.95 < 12 + log(0/H) < 8.2 (Pilyugin 2001). The dif- 
ferent oxygen abundance determination methods and different 
calibration methods may be the cause of the slope of the relation 
slightly deflecting from 1:1. 

The accuracy of the Ar303, Ne302 indicators are compara- 
ble, as can be judged from Fig. [6] The advantage of the Ne3Q2 
method over the Ar303 one is that it does not demand a reli- 
able reddening correction. The most prominent advantage of the 
Ne302 method over other metallicity indicators, such as the R23- 
, N2-, P-, 03N2- methods, lies in its ability to determine the 
oxygen abundance for high redshift galaxies. The Ne302 line 
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Fig. 6. Comparison of the Ne302 and Ar3<93 methods. Ar303 
metallicity is calculated using the calibration of Stasinska 
(2006). Ne302 metallicity is calculated using Eq. (2). Sample 
I is indicated by stars, and Sample II is indicated by crosses. 

ratio can be measurable at the extreme redshift limit of ground- 
based optical surveys. The Ne302, R^, N2, P, 03N2 line ratios 
can be detected within several NIR atmospheric windows for the 
specific redshift z~1.5 (Maier et al. 2006); but when it comes to 
the z~5-6 galaxy, only the Ne302 line ratio can be observed with 
near-infrared instruments. When it comes to the z~ 12 galaxy, the 
Ne302 line ratio can be detected with the James Webb Space 
Telescope (JWST), which uses a Near Infrared Camera and a 
Near Infrared Spectrometer (NIRSpec). 

5. Conclusions 

From the literature and SDSS DR4, we collected a large sample 
of spectroscopic measurements of H n regions and H 11 galaxies 
covering a wide range in metallicity (7.0 < 12 + log(0/H) < 
9.0). The T e values for 210 galaxies (Sample I) in SDSS DR4 
are calculated with direct methods and 2960 galaxies (Sample 
II) with an empirical method. We confirm the existence of the 
correlation between Ne302 index and oxygen abundance and 
have obtained a calibration that can be used as a method for de- 
termining oxygen abundances. The dispersion between oxygen 
abundance and the Ne302 index may come partly from the mod- 
erate depletion of oxygen onto grains. 

Though the Ne302 method is an empirical abundance deter- 
mination method, we believe that also using the Ne302 index as 
a metallicity calibrator presents several advantages: 

1. The Ne302 versus metallicity relation is monotonic. 

2. The relationship between Ne302 and metallicity is not af- 
fected by the chemical evolution effects because the Ne and 
O are produced by the same stars. 

3. The Ne302 line ratio relies on the ratios of emission lines 
that are close in wavelength so do not depend on reddening 
corrections. 

4. This Ne302 line ratio can be measured for galaxies up to 
Z ~ 1.6 using an optical telescope, up to z ~ 5.2 using near- 
infrared instruments on the ground-based facilities, and up 
to z ~ 12 using JWST/NIRSpec; therefore, this flux ratio is 
a promising tool for metallicity studies at high redshift. 
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